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Abstract 



^vg . We study electromagnetic transitions of excited baryons in a deformed oscillator quark 

^^ ■ model, where baryon excited states are described as rotational bands of deformed intrinsic 

\JO • states. We describe all necessary tools to compute transition amplitudes in multipole basis, 

^s ■ which are then related to the commonly used helicity amplitudes. We pay a special attention 

on the sign of the amplitudes as well as their absolute values by computing the photon and 

o> ■ 

Q^ ■ pion couplings simultaneously. We have found that the effect of deformation on the transition 

amplitudes is rather weak. The difficulty in reproducing the empirical amplitude of the Roper 
I ' state is discussed. 
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1 Introduction 

Excited baryons as well as their ground states provide precious information on the dynamics of 
low energy QCD. Not only their masses but also various transitions are particularly important 
for the study of their structure and interactions. Recent and future experiments planned at 
facilities such as TJNAL, COSY, CELCIUS and SPringS are aiming at obtaining more detailed 
information on hadron structure ||l|. 

Based on conventional effective models of baryons, e.g., quark models, one can draw a simple 
but rather unified picture for both ground and excited states 0]. There, constituent quarks are 
placed in a confining potential and are interacting through suitable residual interactions. As 
a consequence, we have an intuitive picture similar to the atomic and nuclear physics; baryon 
states are described as single particle states of three valence quarks. 

An extensive work was performed by Isgur and Karl in the non-relativistic (NR) quark 
model which has been applied to a variety of hadronic phenomena |^], including electromagnetic 
couplings of excited states [^. Recently, the role of the Nambu-Goldstone bosons has been 
emphasized in the context of the chiral quark model |^]. In these models, important physics 
is dictated by valence quarks which are confined inside hadrons and interacting each other 
through residual interactions. In the Isgur-Karl model, one gluon exchanges are introduced, 
while in the chiral quark model fiavor dependent forces generated by Nambu-Goldstone boson 
exchanges are included. These residual interactions yield hyper fine splittings in baryon masses 
with configuration mixing. In general, however, the quality of these models crucially depends 
on the type of residual interactions with considerable amount of parameters. 

In such a situation, we still hope to find a better description for the entire baryonic system. 
In this respect, we have recently found a remarkable systematics in light flavor (n, d and s) 
baryons |^, |^. The systematics applies essentially to all SU(3) baryons including almost 80 % 
of the observed states nominated by the Particle Data Group |^, ^]. The resulting common 
structure of baryon spectra then look very much like a rotational band. This in turn implies 
that excited baryons are likely to be deformed in space. In fact, the idea of deformed baryons is 
not new [|l^, 11|. However, our finding has shown that the picture holds for more baryons than 



it was originally thought. Furthermore, it turns out that a very simple non-relativistic quark 
model with a deformed oscillator potential can explain this aspect at once. This is indeed a 
remarkable fact since the model essentially contains only one parameter. We call the model the 
deformed oscillator quark (DOQ) model as described in detail in Refs. P, M] 

Naturally as a next step, an interesting question arises; how does the DOQ model describe 
various transition processes? This is the issue we would like to address in this paper. For this we 
study the electromagnetic transitions of excited baryons. In principle, we can consider various 
transitions between both excited and ground states. Experimentally, however, only transitions 
from the ground state nucleon to excited states are observed, where experimental amplitudes 



are compiled in the form of the helicity amplitudes A1/2 and ^3/2 |^, [l^i IHi M - 

We compute transition amplitudes first in the multipole basis in analogy with nuclear tran- 
sitions. One significant difference, however, between the baryon and nuclear transitions is that 
in the latter the long wave length approximation can be used as the typical energy scale is a few 



MeV, while in the baryonic case, the same approximation can not be used since the relevant 
transition energy is compatible with the inverse length of a typical baryon size. 

The plan of this paper is as follows. In section 2, we introduce the DOQ model and briefly 
discuss the mass spectrum of SU(3) baryons [^, |^. Then we present wave functions for the 
deformed excited states which are necessary for the computation of transition amplitudes. In 
section 3, we study electromagnetic transitions and derive various matrix elements. We discuss 
carefully theoretical amplitudes that can be compared with empirical amplitudes including their 
signs. Appropriate treatment of the sign requires additional information of the pion coupling, 
which we also discuss. In section 4, we compare theoretical amplitudes with experimental data. 
After discussing the results in the naive DOQ model, we consider two important effects which 
are missing in our previous study. They are diagonalization of the non-orthogonal basis and 
relativistic effects. We pay special attention to the transition of the Roper resonance, since this 
is the channel for which theoretical explanation is difficult. Final section 5 is devoted to the 
summary of the present work. 

2 Deformed Oscillator Quark Model 



In this section we briefly discuss the DOQ model for excited baryons |11, 15|, and present 
necessary ingredients to compute transition amplitudes. 

2.1 Deformed intrinsic states and rotational spectra 

Let us start with the Hamiltonian of the DOQ model: 

^^^"^ = E (1: + t(^>" + ^yV^ + ^"^")) ■ (2-1) 

Here the index i runs over three valence quarks. The masses of constituent quarks m is taken 
about 300 MeV for u, d and s quarks. Practically, in the harmonic oscillator model, only 
relevant parameters are the oscillator parameters uii and the actual value of m is not important. 
Applying the following coordinate transformation 

R = — ^(fl + X2 +X3) , 

v3 

p = ^{xi-X2), (2-2) 

1 

A = —7={xi + X2 -2x3) , 
v6 

the center of mass coordinate can be removed and the resulting intrinsic hamiltonian is given 
as 

i^"" = §1 + ^ HpI + -^IpI + -^IpI) + (p - A) . (2-3) 



The eigenstates of the intrinsic hamiltonian ( p-3| ) are specified by a set of osciUator quantum 
numbers {nx,ny,nz; n^,ny,n^) and the corresponding eigenenergies are given by 

= {N^ + l)uj^ + (A^^ + l)ojy + {N, + l)uj, , (2-4) 

where A^^^ = rix + tl^, Ny = riy + n^, Nz = Uz + n^- Thus the intrinsic energy (|2-4| ) is regarded 
as a function of N^, Ny, Nz and w's. 

For a given (Nx,Ny,Nz), the energy minimization is performed with respect to the defor- 
mation 5iVx, S^y and 6uJz- Unless we have some additional conditions, the variation leads to the 
trivial result; i?™* — > 0, when w^ = ujy = ujz ^ 0. In order to avoid this collapse of the system, 
we impose the volume conservation condition Ux^y^z = uj"^ = const. The validity of such an 
assumption depends on the underlying dynamics of quark confinement. In fact, one may relax 
the volume conservation by, for example, adding a term like B/{uJxU!yUJz)^ to the energy ( |2-4| ) . 
It turns out, however, that the result does not change very much, and therefore we simply adopt 
the condition of volume conservation. 

Now the minimum energy of (2-4) is given by 



^-t = 3{Nx + ly^HNy + I^^HNz + l)'/^a; , (2-5) 

when 

111 . , 

uJx '■ u!y : ujz = -T-: : -T-; : -TT . (2-6) 

y "" Nx + l Ny + 1 Nz + l ^ ^ 

The inverse relation of ( p-6| ) is slightly convenient as it gives the ratio of the lengths of axes of 
the deformed state: 

ax:ay:az = Nx + l:Ny + l:Nz + l, (2-7) 

where ax,y,z are root mean square lengths of the x, y, z directions. The ground state with A^ = 
(Nx = Ny = Nz = 0) yields the spherically symmetric intrinsic state just as in the conventional 
quark model. For excited states iV = 1 ( Nx = Ny = 0, A^^^ = 1 when the z-axis is chosen as a 
symmetry axis), the intrinsic state deforms prolately with the ratio of short to long axes 1 : 2. 
The energy of this deformed state is 3.780lij as compared to 4a; of the spherical quark model. 
For excited states N = 2, there are two cases. One is A^^j^ = Ny = 0, Nz = 2 for the prolate 
deformation with the ratio of short to long axes 1 : 3. The energy is 4.327u; as compared to 5a; 
of the spherical state with (n, /) = (0, 2) or (1,0). Another is Nx = Ny = l,Nz = which yields 
the oblate deformation with the ratio of short to long axes 1 : 2, whose energy is 4.762a;. In 
general, the prolate shape takes the minimum energy for a given A^ = Nx + Ny + Nz, and so, 
in the following discussions we consider only the prolate deformation with the symmetry axis 
chosen along the z direction. We summarize in Table |^ physical quantities of prolately deformed 
states. 

Since the deformed states break rotational symmetry, it should be recovered for eigenstates 
of angular momentum /. This can be performed by the standard cranking method |jl5|. The 
key quantities are the moment of inertia I and the angular momentum fluctuation {P) of the 



Table 1: Various quantities of prolately deformed states. Dimensional quantities are scaled by an 
appropriate power of uj. 
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Table 2: Excitation energies of the N = 2 and A^ = 1 rotational bands. Energy 
values in units of MeV are computed by using u = 644 MeV. M 
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deformed state. In the DOQ model, they can be computed analyticahy. For a prolately deformed 
state of z-axis symmetry (w^ = Wy / w^), they are given by 

^N = ( '-^-^ + '-^^-^ ] = ,,, ;,,i/, [N^ + 2iV + 2) -^ , (2-8) 




(^% = - -l\=N{N + 2), (2-9) 



The rotational band is then constructed for each N: 

Eni = E^§'-^ {l')^ + -^l{l + 1) . (2-10) 

Numerical values for I and (P) are shown in Table [l|, and the rotational energies E^i for the 
N = 2 and A^ = 1 bands are shown in Table y . We emphasize that the energy subtraction 
due to angular momentum fluctuation, the second term of (2-l[l|), is very important to make 



the theoretical masses down close to experimental values. This is particularly important for the 
first 1/2+ excited state (the Roper). 

Coupling the intrinsic spin s of three quarks with the orbital angular momentum I, take 
for example s = 1/2, we consider mass spectra of the N = 2 and A^ = 1 bands for nucleon 
excited sates as shown in Fig. ^. On the left hand side we show the theoretical results: one 
for the positive parity states of I = 0, 2, 4, • • • and the other for the negative parity states 
of / = 1,3, • • • . In the theoretical side, two spin states j = / it 1/2 degenerate when spin- 
orbit coupling is ignored, as experimental data suggest. On the right hand side, experimental 
masses of well observed nucleon excited states with four stars are shown [^ . One exception is 
the 5/2~ state of Di5(2200) with two stars. This state is very likely to form a spin doublet 
with Gi7(2190). We do not list all the states but those which are well identified with the ^8 
representations (s = 1/2) of the spin-flavor group, and are well compared with the DOQ model 
predictions. Our theoretical formula for a fixed N {N = 2 for positive and A^ = 1 for negative 
parity states) corresponds to the rigid rotor approximation. Theoretically, one may consider 
effects from higher rotational bands of A^ > 3 in order to account for the softness of the intrinsic 
state. This is discussed in Ref. |16|. 



The mass formula ( 2-lCl| ) should be compared with excitation energies of baryons. Hence it 



contains essentially one parameter, which is the oscillator parameter uj. It is determined here 
by the average mass splitting between the first excited states of 1/2+ (Roper like states) and 
the corresponding ground states for the flavor SU(3) baryons (for details of SU(3) baryons, see 
the discussion below). The resulting value is a; = 644 MeV Q. Considering the simplicity 
of the DOQ model, it is remarkable that many observed states fit very well to the theoretical 
rotational band. In particular we note that the first 1/2"*" excited state of the Roper resonance 
emerges as the band head of the N = 2 rotational band. 

If the fundamental structure of the excitation spectrum is produced by gluon dynamics of 
quark confinement, we should see a similar pattern in other members of the spin and flavor 
multiplets. In a recent publication, we have analyzed in detail the mass spectrum of SU(3) 
baryons and found that the picture of deformed baryons works extremely well for a wide class 
of baryons |^, ^. We have been able to explain almost 80 % of the observed baryons as they fit 
into the rotational bands of the DOQ model. 



Figure 1: Experimental nucleon mass spectra as compared with the 
DOQ predictions. 

2.2 Three quark states 

Baryon wave functions consist of spatial, spin, isospin and color parts. Here we discuss mainly 
spatial wave functions with definite permutation symmetry and with a good angular momentum. 
Since detailed method is found in Ref. |11], we summarize minimally what we need in the 



computation of transition amplitudes. 

First we consider single particle wave functions in the deformed oscillator potential. They 
are given as products of one-dimensional harmonic oscillator wave functions tpn{x)'- 

i^n^nyn-A^) = ^n, (x)V'nJy)^n, (2^) ■ (2-11) 



Since we are interested in excited states of prolate shape, we set Ux = riy = 0, and write ( p-11 ) 

as '(Aooi.(^) = V'o(a;)^o(?/)V'f (^)- 

The state with a definite permutation symmetry for a three quark system, after removing 
the center of mass coordinates, is given by the product of the single particle wave functions 
V'oon(p) and V'oon'(^) with the constraint n + n' = N. They are either totally symmetric (S) 
or mixed symmetric (MS) f\ . The MS states are further classified into two types of p and A 
symmetries. Here we summarize single particle wave functions and the resulting intrinsic states 
for deformed baryons, ^^'^(p. A), where a represents a permutation symmetry [S, A or p). In 
the following we introduce the parameter a = \fmij and the deformation parameter d = ujx/cuz- 
As we see in the previous section, d becomes A^ + 1 after the energy minimization. However, 
we keep it in the following expressions of wave functions to compute transition amplitudes as 
functions of the deformation parameter d. 



"^For the prolate deformation, there is no totally antisymmetric wave function like vE'^'* = <^ioo(p)0ooi(A)- 
^ooi(p)0ioo(A), since quarks excite along only one direction, where {nx = Uy = 0,nz ^ 0). 



• Single particle wave function 

/ 2\ I 
V'ooo(x) = (^ j ' exp(-^{d^(x2 + y2) + d-h^}) 

/ 2\ I 

V;ooi(f)= f^y y2d-iazexp(-^{d^(x2 + y2) + rf-|^2|) (2-12) 

• Three quark state 

* A^ = (Spherical ground state) : d = 1 

vI/0^(p,A) = V'ooo(/?)^ooo(A) (2-13) 

* iV = 1 (Negative parity excitations) : d = 2 

*i^(p,A) = Vooo(/o)^ooi(A) 

* N = 2 (Positive parity excitations) : d = 3 

*^^(/^>A) = —= (^^oo2(a)V'ooo(A) +V'000(/O)V'002(A) 
^2p(^,A)=V'ooi(A)V'ooi(A) (2-15) 

*^^(P) ^) = ^ (-0002(^)^000 (A) - V'ooo(p)V'oo2(A) 



(2-14) 



Since these deformed intrinsic states ( p-14 ) and ( p-15| ) are not eigenstates of angular mo- 



mentum, we need to project out the states with good angular momentum for physical baryons. 
The projection method is conveniently performed first by expanding the deformed single par- 
ticle state ipoouix) by the wave functions (j)ni,m=o of the spherical three dimensional harmonic 
oscillator: 

V'OOi'(^) = J2 ^nl4>nl,m=Q = ^ C'^l(t>nl ■ (2-16) 

nl nl 

The constraint m = reflects the axial symmetry around the z-axis. Substituting ( |2-16 ) in Eqs. 
( 2-14 ) and ( 2-15| ), we find for three quark states ^: 

^^'^(p,A)= Y. Kln.l.KlM^n,l,{\), (2-17) 

where the coefficients F^'f „ , are given by 

pip _ (^1 fiO 

riplpUxlx ■n-plp'-'nxlx ' 

^UplpHxh ■n-plp'-'nxlx ' 

Tp2S _ 1 f^2 ^0 I r'O r'S ^ 

^riplprixh ~ ^2 V ""'" "^'^ ^^riplp'-'nxlx) ' (2-18) 

^UplpUxh ^uplp'^nxlx ' 

7?2A _ J_ ( (^2 fiQ _ fiQ f<2 

^Uplpuxh /q \ nplp'-^nxh ^Uplp'^nxlx 



The states with definite angular momentum I, m can readily be projected out by operating 
the projection operator [|lj, IS], 



Im 



I d[A]Dl:,oiA)npiA)nxiA) , (2-19) 



where TZp{A) and TZx{A) are the rotation operator acting on the p and A variables, and A 
denotes a set of Euler angles for rotation. Applying Pi^ to the deformed intrinsic state ( ^-17] ), 
we find 

^f^^Pl^^''-=Af Y. Kln,lAlpOlxO\lO)[<t>n,l,x4>n,J"', (2-20) 

where [4>npip x '/'nA^A]''" expresses that (pnpip and (pnxh ^^^ coupled to the state with total an- 
gular momentum {l,m), and (Ip Ix 0\l 0) are the standard Clebsh-Gordan coefficients. The 
normalization constants M in ( p-20| ) are given by 

^''= E I<Wa(^p0^aO|/0)|2. (2-21) 

nplpTixlx 

Here we have performed angular momentum projection after taking variation (PAV) with respect 
to (j's. One could get better wave function through angular momentum projection before 
variation (PBV). We have checked the PBV and found that the difference between the two 
projection schemes is not very significant up to / = 2. 

Finally to close this section, we write spatial-spin- flavor wave functions which are totally 
symmetric. Let us denote spin and isospin wave functions by x and <j), respectively, with a 
superscript p, A or S. For nucleon states (denoted by A^) with isospin 1/2, there are three 
possible states for a given spin and parity j^{P = (—1)'): 

(Ni) \N;[is, i/2y) = ±([^f^,xpyr + [^f',x^y^' 

(N2) \N; [1ms, l/2y) = \ ([^f", x'Y cj^^ + [^f", X^V 'P" 

+ [<i>f\ xy4>' - [^^\ x']^'<^') , 

(N3) \N-[lMS,i/2Y) = ^([^'^'',x'Y<t^' + [^f\x'Y<P^ 

On the left hand side, the notation is such that the orbital wave function is labeled simply 
by the angular momentum I and permutation symmetry a^ lu. The orbital angular momentum 
is then coupled with the three quark spin s = 1/2 or 3/2 to yield the total nucleon spin j. For 
deltas with isospin 3/2, there are two possible states for a given j^: 

(Al) |A;[/s, 3/2P) = [^^■'^^x'n', 

(A2) |A;[/a/5, 1/2F) = -^ ([^f", xT + [^f^', x'^^'^ ' ^^"^^^ 



(2-22) 
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Figure 2: Pion photoproduction which is decomposed into resonance 
and background terms. 



3 Electromagnetic Transitions 

3.1 Helicity amplitudes for the pion photoproduction 

Experimentahy, electromagnetic couphngs of excited baryon are extracted from single pion 
photoproductions (Fig. P). In a resonance dominant model, production amplitudes are assumed 
to be decomposed into the Born terms and resonance contributions. A distinguished feature of 
this process is that not only the magnitude but also the sign of the resonance contributions can 
be determined relative to the sign of the Born terms |12, 13, 14|. More precisely, the sign e in 
the following combination can be determined experimentally: 



{N\H^\N*) ■ {N*\H^\N) = e\{N\H^\N*) ■ {N*\H^\N)\ , 



(3-1) 



where Ht^ and -fT-y are the interaction Hamiltonian for the pion and photon couplings. The 
sign e is then included in the electromagnetic couplings {N*\H^\N). Therefore, in theoretical 
calculations, both the electromagnetic and pion couplings have to be calculated simultaneously. 
Ignorance of the pion coupling part might lead to incorrect results for the sign. 

Usually experimental amplitudes (T-matrix) are presented in the helicity basis, which are 



given as matrices in the space spanned by the initial and final helicity states [12, uM. They are 



also functions of the scattering angles 9 and <f>. The helicity amplitudes are then expanded by 



multipoles as ||l^ 



VAl 



j 



(2j + l)^L4.We 



i(X-fi)ct) 



(3-2) 



Here A = A^ — Aj and /i = Ayr — A/ = — A/ are the initial and final state helicities with A^, Aj, 
Att and Xf being the helicity of the photon, of the initial nucleon, of the pion and of the final 
nucleon, respectively. In the pion photoproduction, A takes four values of —3/2,-1/2,-1-1/2 
and -1-3/2, while fi = —1/2 and -1-1/2. Among the eight components of A^X: four of them are 
independent due to time reversal symmetry. The expansion coefficients A-' , which were called 
the helicity elements in the literature [O] contain dynamical information of the excited states 



10 



of the total spin j. The hehcity elements A^ ■, in its form, however, do not have a definite parity. 



^A 



Parity eigenstates are projected out by the following linear combinations 11^ , 13, !£]: 

^a" = ~7| (^l/2A ~"^-1/2aJ ' (^'^) 

where the total spin j is given by j = In + 1/2, and the parity is P = — (— )^'' for C^ and 

P = (— ) for C^^" . Here !■„■ is the orbital angular momentum of the pion in the final state. 

Resonance parameters are determined by fitting the experimental amplitudes ( p-3D and ( |3-4| ) 

by a simple Breit-Wigner form assuming resonance dominance 



1/2 

Ci-^iW) = - y e ( ^-^ 1 — ^ —^ —r; + background 

^ ^ ■' /-^ \ hn I T/r/2 _ A/f *2 _ „■ MV ^ 



N* 




1^2 _ M*2 _ ^^p 

where T is the total decay width of the resonance A^* , T^ and F^ are the partial width for 
gamma (with the helicity A) and pion emission, k and q are the momentum carried by the 
photon and the pion, respectively. In the second line of (3-5) the CM energy W is set to the 



resonance energy M* , where the single resonance A^* is assumed to dominate the amplitude. 

The relative sign of each term is included in the parameter e = ±1. The helicity amplitudes for 

jP 
the electromagnetic couplings A^ are then defined by the imaginary part of the amplitude: 

< = -iwh)i-,wW'''' ""'''"'"'^ '''''• ^■ = '.±V2.P = -H'-, (3-6) 

where C-^n is the isospin Clebsh-Gordan coefficient for the decay of N* into the relevant vrA^ 
charge state [^]. In this equation, the strength of the pion coupling in C^ is removed by 

— 1/2 

multiplying the factor T^ 

Theoretically, the amplitude C^ (M*) may be computed as 

aW) = ^imH.\N*)^,^-^^L__^M*\H,\N) 

^^^^^ {N\Hn\N*){N*\H^\N)l;. (3-7) 



From this and (3^6), one obtains an alternative expression 



A, = -e (KF) {-^y {N*\H,\N) , (KF) = ^|^ , (3-8) 

where / is an orbital angular momentum of an excited quark in N* , and e is the sign extracted 
from, together with another phase factor (— i) , the pion matrix element as given in Appendix B. 
At the end of this subsection, we summarize selection rules for various amplitudes and excited 
states in Table 0. 



^Note that the expression the resonance propagator (3-5) differs from that given in Refs. O, 041. However, (3-5) 
agrees with, for instance, eq. (2-18) of Ref. Q. 
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Table 3: Relation between various amplitudes and quark model states. 
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Figure 3: A photon coupling to a quark in baryons. 



3.2 Photon couplings 

The electromagnetic coupling is given by the interaction haniiltonian 



e (fxJ ■ A 



-fz-v 



where A is the photon field, and for the current J we adopt the non-relativistic form 



J^^) = [u\{iV -iV)ui + SJ >i{u\aui)\T^ 



(3-9) 



Here we have used overall permutation symmetry of quarks in the baryons, and Ui and Uf are 
the two component Pauli spinors for the initial and final state quarks, respectively (see Fig. ^). 
The isospin matrix r^ is either tq = 1 for the isoscalar current and tj = tj, for the isovector 
current. We will consider relativistic corrections later. 

Helicity amplitudes ^1/2,3/2 are then computed as given in ( p-8| ) : 

1/2 = -e(-i)'(i^F)(iV*(i, l/2)|J.i'+i|iV(l/2, -1/2)), 
I3/2 = -e(-i)'(ifF)(7V*(j, 3/2)|J.l+i|iV(l/2, 1/2)). (3-11) 



A 
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In order to relate the helicity amplitudes with the multipole amplitudes, we consider the mul- 
tipole expansion for the photon field |^ 

Ax = ey^e''^-^ . (3-12) 

Here e\ is the polarization vector of helicity A = ±1. The momentum k is related to the physical 
momentum kp^yg by 



^ ~ \j l^^phys ) ("J'l'^j 



in the convention of coordinate transformation ( |2-2| ). If we choose the z-axis along the direc- 
tion of photon propagation k, the plane wave of a definite helicity ( 3-121) can be expanded by 



multipoles as ||2l[l 



^ = -V2^^i'V2rTl(AiSf^+zl[f) . (3-14) 



nx J ■ 

1=1 

Here the magnetic (M) and electric (E) multipole fields are given by 

A'J\^ = Ukr)YUr) , Y^ = -=i==Lll^(f) , 

A^Hr-) = -^Vx4f(f), (3-15) 

where we have adopted the standard convention for the spherical harmonics |22], and ji{kr) are 
the spherical Bessel functions. 



Substituting ( |3-14 ) in (|3-S| ) and applying the Wigner-Eckart theorem to each multipole 



amplitude, we find the following formulae relating the helicity and multipole amplitudes: 

For j = 1 + 3/2: 



Ai/2 = -e{KF)V2^\l^±^MiMl + l) 



A3/2 = -€{KF)V2^J^^j—-^MiMl + l) (3-16) 



/ + 3 
2/ + 4 

For j = 1 + 1/2: 



Ay, = e{KF)V2^[+^ll±^M{Ml + l)-^^M{El) 



A3/2 = e{KF)V2^[-\I^^M{Ml + l)-^^±^M{El)\ (3-17) 



For j = 1- 1/2: 



Ay, = e{KF)V2^ l+^r—^ M{Ml-l) + \r-^ M{El)] 



21 ' ' V 2/ , 



Ay, = e{KF)V2^l+\r-^M{Ml-l)-Jl-^M{El)] (3-18) 



21 ' ' \l 21 , 



^ The same symbol A is used for the photon field and helicity amplitudes, and e for the photon polarization vector 
and the sign of the amplitudes. Both notations are standard and there will be no confusion in the following discussions. 
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For j = 1- 3/2: 



^1/2 



-e (KF) V2tt 



21 



M{Ml-l) 



A 



3/2 



+e {KF) V27r 



21 



M{Ml-l) 



Writing |A^) -^ \i) and |A^*) — > |/), the reduced matrix elements are defined by 



M{Ml)^i{f\\J-A^ 



TWi 



{f\r^^ 



Ml / ' 



M{El)^i{f\\J.Ar'\\i)-{f\T,\i) 



(E), 



(3-19) 



(3-20) 



Mr/ ' 



where the isospin part is the ordinary matrix elements. In these equations, the phase i is 
introduced to make the matrix elements Ai pure real. Furthermore, in Eqs. ( |3-16 ) - ( 3-19 ), the 
phase factor (— i) in ( |3-11 ) is canceled by another factor i coming from the multipole expansion 

(ED- 



Let us briefly outline how to compute the multipole amplitudes ( 3-2C ). For illustration, we 
consider magnetic transitions Ml. Let us write the reduced matrix element as 

1 



7. aW 



{f\\J-A\ 



2m y^ ^\(^ -^)^i +^ ^(^\^^i)) ■A^'^ 



(3-21) 



where on the right hand side it is understood that the integral and sum are taken over configu- 
ration and spin spaces. The initial and final states are written as Ui and uj in order to indicate 
the operation of derivatives. The isospin part can be treated separately in a trivial manner and 
is not shown here. Now in (3-21) integrating the left derivative (V) terms by parts and using 
the relation V ■ A = and x ■ A^ ' = (the radial vector x appears when the derivative hits the 
spherical initial state which is the ground state nucleon), one finds that the convection (V — V) 
term vanishes. Thus only the spin term survives: 



Mi Ml) = u\aui ■ V X AY"'> . (3-22) 

2m J 

Using the relation between the magnetic and electric fields ( |3-15 ), one arrives at the expression 



^<*"' = 2s; 



I 



21 + 1 



(/||[yZ+lCT]'jZ+l|N) 



+ 



l + l 



{f\m_Mji-i\\^ 



21 + 1 
For electric transitions El, again integrating by parts, one finds 



{f\r,\i) 



M{El) 



1 



liE) 



k 



Uf V Ui- A) -\ Ufdui ■ A 



r(M) 



yji J 2m, -I 

Since the initial state is spherical, the first term can be written as 



Uf V Ui) ■ A 



'(E) 



u\u'i X ■ A 



'(E) 



Then using the following properties of the spherical harmonics, 

jl{kr) 
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kr 



Yin 



(3-23) 



(3-24) 



(3-25) 



(3-26) 



for the first term and 

a ■ 4f ) =a-YL jiikr) = [11, a]'" Mkr) , (3-27) 

for the second term of ( 3-24| ), one arrives at the expression 



M{El) 



v^i(i±iI(/||r,^||0 + i(/|||lVl'i,(.r)||i) 
m kr 2m 



{f\T^\i) , (3-28) 



where \i') = d/dr\i). Equations ( ^-23 ) and ( 3-28)) are the general expressions for electromagnetic 



transitions A^7 -^ N* . Concrete expressions for various matrix elements are given in Appendix 
C. 

4 Results and discussions 

We have computed the matrix elements as given in (|C-13| ) - (|C-22|) for various excited baryon 
states of A^ = 1 (negative parity) and N = 2 (positive parity) band^. Matrix elements are 
then investigated as functions of the deformation parameter d. Note that the actual d's, when 
determined by the energy minimization, are d = 2 for N = \ and (i = 3 for N = 2 bands, 
respectively. We use the oscillator parameter as determined from the mass spectrum: uj = 644 
MeV. For the photon momentum k we use the values computed from the theoretical masses of 
excited states. The use of the experimental values do not change the essential features in the 
following discussions. 

We have computed the amplitudes for all the DOQ model states which contain more states 
than experimentally observed. For a more realistic treatment, we should consider configuration 
mixings due to residual interactions. They could be induced by gluon and meson exchange 
interactions. For instance, a tensor force between quarks strongly mixes the two states of "^Pms 
and '^Pms for 1/2- iV(1535) and Ar(1650) § ||]. 

In Tables ^^, we summarize helicity amplitudes for photon couplings to nucleon and delta 
excitations. The first column is for the present results with suitable deformation (d = 3 for 
positive parity states and d = 2 for negative parity states) and the second column is for the 
spherical limit (d = 1) which corresponds to the conventional non-relativistic (NR) quark model. 
These theoretical values are then compared to experimental data listed in the third column, 
where we indicate observed states which are naively identified with one of the DOQ model 
states. 

As explained in the preceding subsection, we have computed both the photon and pion 
couplings simultaneously. Therefore, there is no ambiguity in relative signs of the photon 
couplings. In the nucleon channels 5/2+ (^Z)^) and 2>/2~(^Pms)i however, the signs of the 
present result do not agree with those of, for instance, the pioneering work by Feynman, Kislinger 
and Ravndal [24|. However, if we included the sign as tabulated by Moorhouse and Parsons [25|, 



^ We do not discuss the transition A(1232) — > N , since in the DOQ model they both appear spherical states and 
the transitions are the same as for the spherical quark model. 
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Table 4: Helicity amplitudes of nucleoli excitations. 









^1/2 






^3/2 






Positive parity 


d^3 


rf=l 


Exp. 


d = 3 


d=l 


Exp. 




1/2+ ^S's 


P 


+115 


+25 


-68 ±5 


- 


- 


- 


Pn(1440) 




n 


-74 


-16 


+39+15 


- 


- 


- 




^Sms 


P 


+ 16 


+17 


+5 ±16 


- 


- 


- 


Pii(1710) 




n 


-5 


-6 


-5 ±23 


- 


- 


- 




^Dms 


P 










- 


- 


- 






n 


+3 


+4 


-5 ±23 


- 


- 


- 




3/2+ ^Ds 


P 


+70 


+ 113 


+52 ± 39 


-22 


-36 


-35 + 24 


Pi3(1720) 




n 


-21 


-34 


-2 ±26 








-43 ± 94 




^Dms 


P 


+57 


+80 




-18 


-25 








n 


-40 


-56 




+ 18 


+25 






^Dms 


P 






















n 


-6 


-8.5 




+ 11 


+15 






Sms 


P 






















n 


-12 


-14 




-20 


-23 






5/2+ 2^5 


P 


+8 


+12 


-17±10 


-44 


-71 


+ 127+12 


i^i5(1680) 




n 


-26 


-42 


+31 + 13 








-30 ±14 




^Dms 


P 


+6 


+9 




-36 


-50 








n 


+ 15 


+21 




+36 


+50 






^Dms 


P 






















n 


-4 


-6 




-17 


-24 












^1/2 






^3/2 






Negative parity 


d=2 


d=l 


Exp. 


d=2 


d=l 


Exp. 




1/2- 'Pms 


P 


+142 


+ 154 


74±11 


- 


- 


- 


5'ii(1535) 




n 


-117 


-126 


-72 ±25 


- 


- 


- 




^Pms 


P 








48 ±16 


- 


- 


- 


5'ii(1650) 




n 


+ 13 


+14 


-17±37 


- 


- 


- 




3/2- ^Pms 


P 


-18 


-19 


-23 ±9 


-127 


-137 


163 ±8 


^13(1520) 




n 


+55 


+59 


-64 ±8 


+127 


+ 137 


-141 ±11 




^Pms 


P 








-22 ±13 








0±19 


£'13(1700) 




n 


-6 


-6 


0±56 


-30 


-33 


-2 ±44 




5/2- ^Pms 


P 








+19+12 








+19+12 


Di5(1675) 




n 


-22 


-28 


-47+23 


-31 


-39 


-69 ±19 
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Table 5: Helicity amplitudes of delta excitations 



Positive parity 


d = 3 


^1/2 

d = l 


Exp. 


d = 3 


^3/2 
d=l 


Exp. 




1/2+ 


Sms 


+ 12 
+11 


+14 
+ 18 


-12 + 30 


— 


— 


— 


P3i(1910) 


3/2+ 


'Dms 

4 c" 


-9 
-12 
-54 


-14 
-24 
-37 


-20 ± 29 


+11 
+20 
-94 


+18 
+41 
-65 


+ 1 + 22 


^33(1920) 

^33(1600) 


5/2+ 


'Dms 


+28 
-8 


+46 
-16 


27+13 


+22 
-33 


+36 
-65 


-47+19 


^35(1905) 


Negative parity 


d = 2 


Al/2 
d=l 


Exp. 


d = 2 


^3/2 
d=l 


Exp. 




1/2- 


Pms 


-35 


-40 


19 ± 16 


- 


- 


- 


^31(1620) 


3/2- 


Pms 


+69 


78 


116 + 17 


+69 


+78 


77 + 28 


^33(1700) 



and an additional phase factor i'^^ for the D-wave component in the multipole expansion (^ 
Q), the present results agree with the results of FKR. The convention of the sign has not 
been considered seriously even in recent publications. In the work of Koniuk and Isgur p|, 
their convention was taken to reproduce empirical one, which has been followed by several 

This situation was partly discussed by Capstick in Ref. ||29||. 



literatures |26, 27, | 

In many channels, the DOQ model (d = 2 or 3) reproduces experimental amplitudes reason- 
ably well including both signs and absolute values nearly to the same extent that the NR quark 
model {d = 1) does. This means that the d dependence is not very strong, which is somewhat 
surprising, since excited states are rather strongly deformed as d = 2 and 3. A possible rea- 
son could be that the number of the relevant degrees of freedom for baryons, which are valence 
quarks, is not very large. If we, however, look closely at the numbers, the positive parity baryons 
have slightly stronger d dependence than the negative parity baryons, as d is larger for the for- 
mer. Typically as it is seen in the channel 3/2+ ^Ds, AP'°*°°(d = 1) = 113 x IQ-^ GeV-^/^ ^nd 
^1/2 °°(^ = 3) = 70 X IQ-'^ GeV~^/^ as compared with the experimental value Ay°^ °'^(exp) = 
52 + 39 X 10-3 GeV-i/2^ 

Theoretical predictions differ most significantly for the first 1/2+ excited state for the Roper 
resonance; not only the magnitude but also the sign do not agree with data. The relevant 
amplitude is ^1/2 or Ml. Indeed the discrepancy in the sign is puzzling, if the 1/2+ excitation 
is simply considered as a radial excitation of the ground state nucleon. It does not matter 
whether it is a single particle or collective excitation. Due to an S-wave nature of the orbital 
motion, spin structure and hence the magnetic moment of the excited state should be the same 
as the ground state. Now we know that the magnetic moments of the nucleon, the delta pO] and 
the Ml transition of the delta [^ are reasonably well explained by the quark model. Therefore, 
one would naively expect that the magnetic transition of the 1/2+ excited state were explained 
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as well as these transition amplitudes in the above. 

The discrepancy in the sign is known for a long time and various ideas to overcome the 



problem have been proposed. In Ref. |26], deformation effects was considered to be one candidate 
of such. However, their conclusion was based on the wrong sign of the amplitude. In fact, the 
deformation is to increase the amplitude but with the wrong sign remained. Since the problem 
in the Roper resonance is in itself an interesting question, we discuss two important effects in 
the next two subsections. 

In the end of this subsection, we comment on the sign of the amplitudes for the A^(1440) 
and A (1600). The latter is interpreted as the Roper resonance of the delta. Let us look at the 



formulae for the nucleon ( C-13 ) and the delta for d = 1 when the sum in the radial integral 



reduces just to a single term. We have 



Ml{N;^S's) = ^Ml{A;^S's) 



3 

472^ 
k 1 



2^ ^^Fofio(00|jo|10) • (isospin part) . (4-1) 

From this, after taking into account the sign of the pion matrix element (see Table g in 
Appendix B), we observe that the signs of the transitions for A^(1440) and A(1600) are different. 
This fact is correctly shown in Tables H and 0. In literatures, however, the signs of the two matrix 
elements are found to be the same [29, |3^ . 

4.1 Diagonalization of the 1/2+ states 

In the DOQ model, the two 1/2^ states for the nucleon and the Roper, '^Sg{N = 2) and 
'^Ss{N = 0) are not orthogonal to each other, since they are eigenstates of different Hamiltonians 
with different deformation: d = 1 and 3. Generally, those states with the same spin parity j but 
with different N are not orthogonal. Explicitly, denoting \^S's{N = 2)), \'^Ss{N = 0)) = |2), |0), 
they have non-zero overlap 

im-Foioo^^- (4-2) 

In other words, the two states |0) and |2) form a non-orthogonal basis. Hence we need to 
reconstruct orthonormal basis. 

Using the non-orthogonal states |0) and |2), the eigenvalue equation takes the following 2x2 
matrix form: 



{H-EN)'^ = 0, (4-3) 



H= !° '!" h N=r :]^ ^ = \ :;:: i-i : i- (4-4) 



where 



£Eq2 E2 J \ e 1 

Here the diagonal components of the energy matrix are given by ( 2-10]) : 





E^ = {N\H\N) = ET - :^{1')n , (4-5) 
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while the off diagonal components are, for simplicity, taken as the left-right average: £'02 
{0\H\2) = {Eq + £'2)72 • Eigenvalues and the corresponding eigenstates are given by 



En,r = EoTil- 1)A , (4-6) 



and 



\N) = ^ (g^|o) + (i-^)|2)), 

V2(i - Vi^^) 

\R) = "^ (-g^|o) + (i+^)|2)). (4-7) 

^2(1 + Vl^^) 

where 7 = l/\/l — e^, A = (£'2 — Eq)/2. Numerical values for eigenvalues in units of lo are 

Eo = 3.0 , E2 = 3.751 , e = -0.538 , (4-8) 

and therefore, we find 

En = 2.930 , Er = 3.821 . (4-9) 

As expected, the mass of the nucleon is lowered and the mass of the Roper is pushed up, where 
the order of the effect can be estimated as 

{Er - En) - {E2 - Eq) ^ 0140 ^ ^ ^^^ 
E2-E0 0.891 

We find that the effect on the mass difference is of order e^. In contrast, the wave functions 
(orthonormal states) are modified substantially: 

|iV) = -hl.l39|0) -F0.332|2), 

\R) = -h0.332|0) + 1.139|2) . (4-10) 



Using the wave functions (4-7), transition amplitudes are modified as 



{R\M\N) = -sgn{e)r^{0\M\2) + ^^Y^) (^^'-^l^^ + ^2|-^|2)) • (4-11) 

Before the diagonalization only the first term with the unit coefficient (e — > —0) contribute. 
The new terms, the second and the third terms, have the opposite sign to the first term, and 
therefore, the effect of the diagonalization is to reduce the amplitude. Substituting numerical 
numbers for the e and matrix elements, we obtain 

^^/2 = 153 - 60 - 57 = 36 GeV"^/2 , 

AI/2 = -102-h40-h38 = -24GeV"^/^ (4-12) 

We see that the diagonalization affects more on transition amplitudes than on masses. This is 
not surprising, since the effect is of order e^ for masses, while it is of order of e^ for transitions. 
We note that the modified values ( [4-12 ) are now rather close to those oi d = 1. 
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4.2 Relativist ic effects 

The importance of relativistic effects was first empliasized by Kubota and Ohta already twenty 



years ago |33|. Further extensive studies were also performed by Capstick and Keister |28, 34 1. 
It should be particularly so for transitions where the leading order contributions in the long 
wave length limit are forbidden due to some trivial selection rules. The Ml transition to the 
Roper resonance is a typical example of such; the Ml operator with a spin flip leaves the spatial 
wave function unaffected and therefore the matrix element vanishes due to the orthogonality 
between the initial and final spatial wave functions in the limit of zero momentum transfer, 
A; — > 0. Hence the non-zero matrix element at finite k is of order 0{k) where the relativistic 
corrections start to occur. 

So far the fully relativistic treatment has not been achieved yet, since it inevitably requires 
field theoretical method for interacting light quarks. In literatures, instead, a classical method 



of 1/m expansion has been often considered p9| , P2| , 33, ^, ^]. We should note that this 
method suffers from a fundamental question whether 1/m expansion converges or not, when 
m ~ 300 MeV and q ~ 500 MeV. Although there is such a difficulty, still we think that the 
leading correction provides us with some fiavor how relativistic effects appear. 
The leading order contribution of the electromagnetic hamiltonian 

H^ = -e Sx-ip^Tp-A , (4-13) 



leads to the standard form of (p-10). The next to leading order term is of our concern here. 



There are two terms which are so called spin-orbit and non-additive terms. Here we consider the 
spin-orbit contributions in order to see the importance of relativistic corrections. Performing 
the Foldy-Wouthuysen transformation up to order 1/m? , we find the relevant terms: 

H^ ~ ^u\ (~2A -p- ia -k X a) Ui-\ :ru\ (a-kxA-2a-Axp)ui. (4-14) 



2m"/ V~"^ /^-- --^J^^^sm^'^f 
Here k is the photon momentum and p the momentum of the initial nucleon (acting on Ui). In 



deriving ( 4-14 ), we have used the condition for the real photon k ■ A = 0. 

These relativistic corrections are computed using m ~ 300 MeV in the DOQ model after the 
diagonalization is performed. We find, after including the corrections, the total amplitudes: 

Al^2 = -13 X 10^3 GeV-^/2 ^ 

A1/2 = 9 X 10^^ GeV-^/2 . (4-15) 

As in previous works, we find that also in the DOQ model the relativistic correction changes 
the sign of the amplitudes. Although the DOQ values in (4-15) are larger than those of the 



traditional quark model, they are still too small as compared with experimental values. We 
note the similarity in the behaviors of the numbers in the conventional and the DOQ models 
when relativistic effects are included. In both cases, relativistic corrections are very large. If 
so, there is an essential question on the validity of the 1/m expansion. This is a very important 
problem for the non-relativistic quark model, although the model seems to work well for many 
phenomenological aspects. . 
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4.3 Limits in a non-relativist ic treatment 

Here, focusing on the electromagnetic transition of the Roper, we demonstrate that within a 
non-relativistic treatment of a naive quark model it is not possible to solve the sign problem of 
the amplitude. Although the statement is rather trivial, it is useful to point it out here. 

As emphasized in the preceding discussions, we need to compute both the electromagnetic 
and the strong (pion) couplings when discussing the sign. To leading order in the non-relativistic 
expansion, transition operators for these two couplings are given by ( 3-10| ) and ( [B-^ ) . Let us 



assume, as in the naive quark model, that both the nucleon and the Roper is dominated by the 
orbital wave function of / = 0: |A^) ~ |Roper = i?) '^ |[0, 1/2]"*^'^). Then, it is not difficult to 
show that both the couplings contain the spin matrix element {a) : 

{N\H^\R) • {R\H,\N) ~ {N\\a\\R) ■ {R\\a\\N) ~ \{N\\a\\R)f . 

Therefore, it is apparent that within this treatment the sign of the amplitude can not be changed. 
We need more delicate treatment for a better understanding of transition amplitudes. 

5 Summary 

In this paper we have studied electromagnetic transitions of excited baryons using a non- 
relativistic quark model with a possibility for excited baryon being deformed. The main purpose 
of the present work is to test the success of the deformed oscillator quark (DOQ) model for the 
masses of flavor SU{3) baryons by examining electromagnetic transitions. Experimentally, tran- 
sitions from the ground state have been observed in pion photoproductions. 

We have derived all necessary formulae for multipole amplitudes, which are transformed 
into the conventional helicity amplitudes. In the comparison of the theoretical amplitudes with 
experimental data, we have paid a special attention to relative signs among amplitudes by 
computing the pion couplings explicitly. 

Our main interest is how electromagnetic amplitudes are influenced by the spatial deforma- 
tion. We have seen, however, rather small dependence on deformation, which is typically less 
than 50 %. To this order, there are many other effects which modifies theoretical predictions. 
Furthermore, experimental amplitudes usually contain even larger ambiguities due to difficul- 
ties in the analysis. Other types of transitions might be useful to investigate. They could be 
transitions between excited states through photon or pion emission. 

Another interest has been paid to the transition to the Roper resonance. The discrepancy 
between theoretical predictions and experimental data both in magnitude and sign should be 
considered seriously. It does not seem that there is a simple explanation so far, if we recall 
the success of the magnetic moments for the ground states of the nucleon and delta. Further 
investigations for the Roper resonance will be an interesting subject. 



Appendix 
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A Helicity formalism 

Empirical amplitudes for the pion photoproduction is presented by the helicity coefficients C, 



of (S^S) and (324). In order to compute them we explore basic ingredients for the helicity 
formalism, following the classic paper by Jacob and Wick ||l^ . Let us assume that the relative 
momentum between the initial photon and the nucleon is along the z-axis with their helicities 
being A-^ and Aj. In the center of mass system, the total helicity is given by A = A^ — A,. This 
is equivalent to the z component of the angular momentum in the spherical basis. Let us write 
the initial state as 

|A;,6i0,A) = |/c,00,A), (A-1) 

where Q^ = 00 are the polar angles of the momentum k. Similarly we write the final state for 
the pion and nucleon system as 

\q,e(t>,fi), (A-2) 

where 9(1) are the scattering angles of the pion, and fi is the helicity of the final state. 
The helicity states ( A-l| ) and ( [A- 2] ) can be related to the spherical states through 



\p, 94), \) = Y1 V ^^^™a(<^' ^' -0)|P, i"^, A) , (A-3) 

jm 



and the inverse relation 



\p,jm,X) = ^1^^ J dnDi^^icl>,9,-<P)\p,9<p,X) 



(A-4) 



In these equations, Dmn{ci, (3,^) are the Wigner's D-function |22]. For rotation of a vector, the 



angle 7 is redundant and can be fixed as a = —7 = (j). The integral measure is then given by 
dQ = d{cos 9)d(j), and so J dQ = Att. 

Applying ( [A-3| ) to ( [A-lD and ( A-2| ), we can rewrite the helicity amplitude (3-2) 



A^x{94>) = {q,9c|),^l\A\k,m,x:, 



= Y.\j'^^'^'^-^^r:^.{(t>.e,-<t)){q,jm,i,\A\k,j'\\). (A-5) 

jj'm 

On account of the conservation of angular momentum, j' = j and m = A, we obtain 

{qM.^^\A\Km,\) = ^Y.^2j + l)D{*^{ct),9,-<|,){q,j\^l\A\Kj\\) , (A-6) 

i 

which corresponds to ( ^-2|) . 

It is instructive to re-express these relations using spherical harmonics. For this we consider 

|fc,00,A = l/2) ~ e*-^xi/2 = 47r^i'-ll,„j£)y,:^JA:)x(A:r)xi/2, (A-7) 

l-rrTfl-TT 
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where Xi/2 is the spin up state for the nucleon. Using Y^ ^^{k = z) = y^(2/^ + l)/47r 5m^o aiid 
recouphng the angular momentum, we find 

|A:,00,l/2) ~ V2^ j;v/2jTT(.^--i/^3^;7;f -.^■+^/^3^'t/f ) , (A-8) 

i 

where 3^,'^ is the spinor harmonics including the spherical Bessel functions ji^{kr): 



jm 



y\^m = ^{km-^i 1/2 fi\j m)^i^^^^ Xi. , (A-9) 

where ^im = yimjiif^i^)- Similarly, we find 

\k, 00, -1/2) ^V2^Y1 a/^JTi {i^-'/'y^^li' + ^^'/'yjtT) ■ (a-io) 

j 

The states with a definite angular momentum can be obtained by applying the projection 
operator: 



\k,jm, ±1/2) ~ ^'^^ J dnD'^^^^^{cP,e,-cl))R{<P,9,-cP)\k,00,±l/2) 



A-K 



■^-'/'y^;li'±i^^'/'yitT)- (A-ii) 



~ ^/2 V J' 1/2 ^'' -^ii/a 

Here R{4>, 0, —cj)) = exp(— iaJ^) exp{—if3Jy) ex.p{—i^Jz) is the rotation operator, and in deriving 
the second equation we have used the definition of the D-functions. Obviously, in these states 



( A-11 ) parities are mixed. States with a definite parity can be obtained by 



\p,jm,±) = ^{\pjm, +1/2) ±\p,jm, -1/2)) =±A7,P^^/^yff/\ (A-12) 

These are the eigenstates of parity P = !]{—)'"' = T]{—y^^''^, where rj is the intrinsic parity of 
the system. 



B The pion matrix element and the sign e 

Let us start with the relation ( |3-6| ) : 

= (•■■){g,jX.±\H,\N'(jX))Gf (N'{jX}\H^\k,j\,\)C^„. (B-2) 

We need to compute the pion matrix element {q,jX,±\H.,r\N*{jX)) in order to account for the 
sign of the whole amplitude. 

To do so, it is convenient to have an explicit form for the final state as given in ( [A-12 ): 

y 

J (B-3) 

± i-iy^ Y.(ln -m 1/2 X + m \ j A) (-)™ {xx+mm^m . 



{q,jX,±\=(±i'-\y^l)y = Ui'^Y^il^m 1/2 X\-mj)^i 

\ m 
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In this equation, l^^ = J T 1/2 is the orbital angular momentum for the pion relative to the 
nucleon. The interaction hamiltonian for the pion-quark coupling is given by 



Ht^ 



2m 



(B-4) 



where the derivative Vyr operates to the pion wave function. Sandwiching this hamiltonian with 
(|B-3| ) and the resonance state |A^*(jA)), we find after some algebra 

{q,jX,±\H^\N*{jX)) 



-^4 9Q -z, 
2m 



y^(l/2|l[«,.«.|'-|l.> + y3Z{l/2||[S,.-,.l'-|U) . (B-5) 



In fact, one of the two terms is relevant depending on the orbital angular momentum Iq = I of 
an excited quark in N*. Due to the pseudoscalar nature of the pion-quark coupling, it follows 
that I = Ijr ±1. When j = I — 1/2 or j = / — 3/2, the first term of (|B-5|) survives, while when 
j = 1 + 3/2 or j = / + 1/2, the second term does. 

For later convenience when performing actual computations, we summarize relevant pion 
matrix elements according io j = I + 3/2, / + 1/2, / — 1/2, I — 3/2: 



1. When j = 1 + 3/2, / + 1/2 



It,=1+1 



J = In ±1/2 , 



{q,jX±\H^\N*{jX)) = ^4vrgi'+yl±l(l/2||[y,a]'+i| 



(B-6) 



2. When j = 1- 1/2, / - 3/2 



u=i-i 



j = ln± 1/2 , 



(,,,A±|F.|Ar*(,A))=±4vr^.'-^ ^^_^ 



/ 



(i/2ii[yzcT: 



\i-i\ 



Writing the pion matrix element as 

{q,jX ± \H^\N*{jX)) ^ -i{-iy e \{q,jX ± \H^\N*{jX))\ , 



(B-7) 



(B-8) 



the sign e is given as summarized in Table ^. In the table, t] is the sign of the radial integral 
^^ F(j^^;(00|jz|n/) . The notations of this equation is defined in ( C-4 ). 



Table 6: The sign e extracted from the pion matrix elements. 





Nucleon 








Delta 






j 


e 






j 


e 


5=1/2 


1 + 1/2 


-flS or A 


S = 


= 1/2 


1 + 1/2 


+rix 




1-1/2 


+'riS or A 






l-\/2 


~V\ 


5 = 3/2 


1 + 3/2 


+m 


S = 


= 3/2 


1 + 3/2 


+Vs 




1 + 1/2 


-m 






1 + 1/2 


-Vs 




1-1/2 


-v\ 






1-1/2 


-Vs 




1-3/2 


+rix 






1-3/2 


+VS 



ria ^sgn [j:^F,%-^mMnl)] 
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C Computation of multipole amplitudes 



Using the wave functions for N* as given in ( 2-22| ) - (2-23) one can compute multipole amplitudes 



in a straightforward manner. For a given excited baryon, there are three possible amplitudes: 
M.{M l + l), Ai{M I — 1) and M.{El). To illustrate how actual computations can be performed, 
we consider one example for M{M l + \) for the transition |A/^(ground state)) — > |A^*; \ls-, 1/2]-')- 
In the following expressions, the factor 3 of (3-10) is not included, although it is in our 
numerical results. Replacing I ^ I + 1 vn ( 3-23| ), only the second term survives due to the 



matching of orbital angular momentum. Since the initial state for the ground state nucleon is 
spherical, / = 0, the orbital angular momentum of the absorbed photon and that of the excited 
baryon must be the same. Thus we have 

M{Ml + l) = -Ay|±|(^([.i,f^VF||[y^^]'||.I.of^^P)^^(p) 

+ ([^rVFII[y^^]'ll^?lox')r,(A)) . (C-1) 

Here we have included the isospin matrix elements which are defined by 

r^(p) = mr^\r) , r^(A) = (0^|r^|0^) , (C-2) 

and the radial integral Jr'^dr is implicit. Now the reduced matrix elements can be computed 
in a straightforward manner using the decomposition theorems, and we find the result 



^ ^ 4mV2/ + 3-^ U + l 1 1/2 



47r 



^Y.^nm {nl\jim V6 (r^(p) - \r,{\)) , (C-3) 



where / = \/2/ -|- 1. Here we have introduced the notation for the radial matrix elements: 

{nl\ji\m) = f r^drRnijiRoo , (C-4) 



where Rni are the radial functions of the harmonic oscillator. 

All other matrix elements can be computed in similar ways. Here we summarize the results 
for all amplitudes. Let us introduce the following notation for the sum of the radial matrix 
elements: 

P,^- = -^Y.^oTniinl\jim, (C-5) 



^^' = -^Y.^^Oniinl'\jim, (C-6) 



and the combinations of isospin matrix elements 



r±i/3 = niip) ± ^M^) ■ 
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We checked the behavior of the coefficients F numerically and found that they converge to zero 
quickly as n is increased. In the present work, we take the sum over nine terms in ( |C-5D and 
(|C-6| ) to obtain the accuracy of more than five digits. 

The matrix elements for various excited states N* are summarized as follows: 

Nl iV* = |[/s, l/2p): 



MiMi + i) = -^J^jm\y,l { /,JxP,^^^r_,/3 (C-7) 



4mV2Z + 3 1/ + 1 1 1/2 



Mm = -^^^(-1)' J QP T,, + ^ri l^f { J^l P,-- V6 T_V3 (C-9) 

N2 N* = \[ImsA/2V)- 

This case is obtained simply by performing the following replacements in the results of 
(Nl). 

• Multiply l/\/2 as an overall factor. 

• Change the sign of one of the isospin matrices: +r^(A) -^ — t^(A). 

• Replace the superscript of the coefficient F: F^^^ —>■ F^^^ and hence Pi^^ -^ pA^A 
and Q^^ ^ Q^\ 

N3 N* = \[lMS,m')- 



M{M i+i) = -A /ill j rn { ^/^ ^ ^Ip,^^a^(a) (c-10) 

^ ' 4mV2/ + 3-^ h + 1 I l/2j ' V3 '^^ ^ ^ ' 



>1(M/-1) = A /_LA(_l)J+3/2 3^-^l 3/2 J MpiVA4 

^ ^ 4mV2/-r ^ ■' \l-l 1 l/2j ' ^/3 '^^ ^ ^ ' 



M{El) = -^jir'^ ^ ^ U^^A 4 
^ ^ 4m^ [ / 1 1/2J ' V3 

The decays of delta excited states Al and A2 are obtained by the following manipulations. 

Al A* = |[/a/5, 1/2]^): In the results of Nl 

• Pick up the r^(A)-term and replace r^ -^ {(J)^\t^\(J)^) . 

. Replace F^,^, by F^i- 
A2 A* = \[ls, 3/2]^): In the results of N3 



• Replace r^ -^ ((/)^|r^|(/)'^) 
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• Replace F^^f, by F^^^, . 

• Multiply \pl as an overall factor. 

For practical purposes, it is convenient to write the spin of excited baryons j by / (j = /it 1/2 
or / lb 3/2), so that the 6-j symbols are computed explicitly. Results are: 

Nl.a iV* = |[/5,, l/2]'+i/2), allowed transitions: M / + 1 and F /, 

M{M / + 1) = - Ay^(n:^ pfS j._^^^ (C-13) 

S rp , k^/2l J^g 

2m -' ^+^ + ^^^^ '-"'' 



M{E I) = ^^^^ Qr T,r + ^ P.-^ T_,/3 (C-14) 



Nl.b N* = \[ls, l/2]'-i/2), allowed transitions: M / - 1 and F /, 

k 



M{Ml-l) = —72(1^ P,^* r_i/3 (C-15) 



4m 



A^(F = - VW^)l QMS T^^ + k^MTT) ^^s T (C-16) 

Im 4m ' 

N2 The results for this case, where N* = \[ls, 1/2]'^^'^) are obtained by making the replace- 
ments as described before. 



N3.a N* = \[Ims, 3/2]'+^/^), allowed transitions: M / + 1, 

N3.b N* = \[Ims, 3/2]'+i/2), allowed transitions: M / + 1 and P /, 



^('^' + ^' = -3^\/|l^^"^^''« <«■'«' 



A<(£0 = -^\/^Frr„(A) (C-19) 

N3.C N* = \\Ims, 3/2]'-i/2), allowed transitions: M / - 1 and F /, 



k_ {l-l){l + l) ^^, 
3m Y 2(2/ -1) 



M{M / - 1) = -£:\l ' ^,;;^ 7,^ Pf"^ T^(A) (C-20) 



A^(^0 = ^V\^^''''"''(^^ (C-21) 



N3.d iV* = |[/m5, 3/2]'-3/2>)^ allowed transitions: M l-l 

k [2 1-1 



^(*"-" = 2^Vt^lr=T^"^''<'' '^-''' 



Similarly, we can derive the formulae for deltas. 
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